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Abstract 

Global harmonic coordinates for the complete Schwarzschild metric are found 
for a more general case than that addressed in a previous work by Quan [4]. The 
supplementary constant that appears, in addition to the mass, is related to the stress 
quadrupole moment of a singular energy-momentum tensor. Similar calculations are 
also carried out in q-harmonic coordinates. 

1. - Introduction 

Since its advent, the theory of General Relativity has been used broadly to handle 
harmonic coordinates for several purposes. In many cases, such coordinates are used for 
reasons of mathematical comfort [1], while other times they have been used with the 
purpose of attributing them a special meaning [2]. In particular, they are considered 
to be a natural generalization of Euclidean Cartesian coordinates. Nevertheless, before 
recently no global System of harmonic and asymptotically Cartesian coordinates has been 
offered for a stellar-like model; that is, a System that is well defined everywhere (inferior 
and exterior) and with a C^ metric on the surface of the star (admissible coordinates 
in the sense of Lichnerowicz [3]). One example of this has been provided by Quan Hui 
Liu [4]: the complete Schwarzschild metric (inferior and exterior) with /iq = 8/9 (/Uq 
being the dimensionless quotient between two times the mass of the star and its radius) . 
In this case, it should be noticed that this special value of the fiQ parameter leads to a 
simplification of the calculations but, unfortunately, it also provides a model of a star 
with divergences in the pressure at its center. This programme forces one to introduce 
two new constants (in addition to the mass) into the metric: one of them, Qext , at the 
exterior and another, Qint , in the inferior. These constants, which are defined with a 
length dimension in the above paper, prove to be proportional to the radius of the star, 
and according to the author are physically meaningless. 

The aim of the present article is dual: on one hand, we shall develop Quan Hui Liu's 
Programme but in a more realistic and general case for the ßo parameter (hq < 8/9). The 
Solution of the problem for an arbitrary value of this parameter leads to a not well-known 
Heun differential equation [5] instead of the simple hypergeometric equation that appears 
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in the limiting case fxo = 8/9. On the other hand, by use of the linear approximation of 
vacuum and spherically Symmetrie Einstein equations we show that the exterior constant 
Qext is closely related to the stress quadrupole moment of the source. 

In Section 2 the complete Schwarzschild metric is written in Standard Droste 
coordinates simply to set the notation to be used, and the Darmois matching conditions 
[6] are briefly reviewed by writing the first and the second fundamental forms relative to 
the surface of the star. 

In Section 3 the change of coordinates is found, from Standard ones to general 
asymptotically Cartesian harmonic coordinates, two conditions being fulfilled: the new 
System of coordinates should preserve the spherical symmetry and diagonal structures of 
the metric (written in associated polar coordinates) and, in addition, the function relating 
both sets of coordinates must be on the surface of the star. These conditions allow 
US to express Quan's exterior and interior constants in terms of the change functions 
evaluated at the boundary. We then prove that the new coordinates are admissible 
coordinates in the sense of Lichnerowicz [3], i.e. the metric is on the surface. 

Section 4 is devoted to interpreting the exterior constant, Qext- To do so, we first 
write the multipole expansion of the exterior Schwarzschild metric in the new set of 
harmonic coordinates found earlier (using the inverse of the radial coordinate, up to Order 
five, as the parameter of the series). That expansion has already been shown, although 
in a slightly different way, by one of us [7]. Secondly, the more general Solution of the 
vacuum Einstein equations with spherical symmetry in the linear approximation is made 
explicit and hence a suitable comparison with the previous expansion can be established, 
providing arguments for understanding the significance of the new constants appearing in 
the metric. Finally, it is checked that this Solution also comes from some Singular energy 
momentum tensor whose stress quadrupole moment proves to be "proportional" to the 
exterior constant of Quan (which means that Quan's constant is the factor multiplying 
the tensorial expression of that multipole moment). 

In Section 5 we briefly discuss the results obtained for the same topic by using 
the so-called q-harmonic coordinates, a variety of harmonic coordinates introduced by 
L. Bei [8]. Most of the contents included here have already been obtained by J.M. 
Aguirregabin'a [9] for the exterior case and by P. Tcyssandicr [10] for both the exterior 
and interior cases. Those results, as well as somc new other ones, are introduced to 
make a comparison with the harmonic scenario and because until now they have not 
been published in their cntircty. 

The paper is completed with an Appendix which is devoted to explaining the Heun 
differential equation; this has been included because of the rather odd and at the same 
time fundamental equation involved. We hope that with this the reader will find it easy 
to understand Section 3. 
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2. - Complete model of the Schwarzschild metric 



The interior Schwarzschild metric has the foUowing expression in Standard polar 
coordinates {t ,r ,9 



dsj = - 



dt'^ + (r) dr^ + r^{de^ + sin'^e d<p'') (1) 



where the radial coordinate r is assumed to be restricted to the interval [0,ro], tq being 
the radius of the star. The foUowing notation is used: 

11 

lir) = 1 - , j^ = -XP iro<L) , x = Stt (2) 

p being the energy density (constant) of the star and where we arc dcaling with 
geometrized units {G = c = 1). As is known, the pressure of the model is given as 
a function of the radius by the foUowing expression 

pw=. n';'-^*y, (3) 

372 (ro) - 72 (r) 

which provides the foUowing values for the pressure at the surface and the center of the 
star respectively 

Piro) = 

p{0)=Pc = P 1 

372(^0) - 1 

Since the value Pc of the pressure at the center of the star must be finite, it relates the 
radius of the star tq and the density parameter L by means of the foUowing restriction 

37^(ro)-l>0 ^ (5) 

The exterior Schwarzschild metric is written in the same set of Standard polar 
coordinates as foUows (r > tq) 



dsl = - (^1 - + (^1 - + r\de^ + sin^ d^^) (6) 

where m represents the mass of the star, which is related to the density p and the radius 
ro by the known expression 

43 / 2m 8\ _ 



The above expressions for the Schwarzschild metric (interior and exterior) are such 
that the Darmois matching conditions [6] are automatically fulfiUed on the surface of the 
Star (E : r = ro). These conditions -that is, the continuity of both the first and the 
second fundamental forms- will be symbolically denoted as foUows 

I: Me = 

, {a,b,... = 0,2,3 = t,e,<p) (8) 

II: [Kab]^ = 

In what foUows we make explicit the expression of the fundamental forms to show the 
requirements provided by (8). 

• The Grst fundamental form is 

hab{p) = 9aß[x{p)]e2/e^/ (9) 
e^^ being the tangent vectors to the surface S 

e-/(/) ^ ^ (E: x'^p'^t, x'=p' = e, x'^p' = ^) (10a) 



which leads to 





t , 


x^ = 


9tt 

























(106) 



i.e., Qu must be continuous, grr being unrestricted (note that the coordinates satisfy the 
Euclidean sphere condition everywhere: g^^ = geesin^O = r^sin^ö). We shall see below 
that the continuity of g^r is imposed by the continuity of the second fundamental form. 

• The second fundamental form is 

^ + rVa/</j = -ra6 (11) 

where the notation Zq, = da{r — tq) has been used to denote the normal vector to the 
surface S 

/drgtt \ 
{Kab) = ^ö"' drgee (12) 
\ drg^^J 

It may be deduced from (12) that both g^r and drgtt must be continuous. As a summary, 
only the discontinuity of drgrr is allowed, which foUowing (1) and (6) turns out to be 
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[drgrr\Y. = drOfri'^o) - OrOrriro) = -^J2^ ^('^o) (13) 

It is then obvious that the set of Standard polar coordinates {t , r , ö , describing the 
complete Schwarzschild model are not admissible in the sense of Lichnerowicz, since the 
metric is not C^. 



3. - Global asymptotically Cartesian harmonic coordi- 
nates 

The aim of this section is to find a global System of harmonic coordinates for the 
complete Schwarzschild model, which must be asymptotically Cartesian and such that 
the components of the metric are C-^ at the boundary E : r = tq. Let us consider the 
generic form of the Schwarzschild metric in Standard coordinates (1),(6) 

ds^ = gu{r)df + grr{r)dr^ + r^{dd^ + sui^ddip^) (14) 
We now perform a change of coordinates 



(16) 



{t, r, e, ^} = {x"} {i, X, y, ~z] = (15) 
satisfying the foUowing conditions 

Jim g^g = r]aß = diag(-l, +1, +1, +1) 

which represent the harmonicity condition and asymptotic cartesian behaviour respec- 
tively, and where f = ^/ x'^ + y'^ + z'^ . 

The easiest way to solve the problem is by using polar coordinates {t, f, 6, (p} 
associated with the harmonic ones {x'^} 

X + iy = f sin e'^'^ 

(17) 

z — r cos 9 

We see that in order to preserve the diagonal structure as well as the structure of the 
spherical symmetry, the change of coordinates (15) should be as foUows 

i=t , f = /(r) , 9 = 9 , ip = ^ 

lim = 1^ (18) 

r^oo r J 
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Henceforth, equation (16) turns out to be the foUowing second-order linear difFer- 
ential equation for the function /(r): 



Let US now see the form of this equation in the interior and at the exterior of the star. 

A) The exterior problem 

For the expressions of the exterior Schwarzschild metric given (6), equation (19) is 
written as foUows 

r (r - 2m)/"(r) + 2(r - m)/'(r) - 2/(r) = (20) 
whose general Solution is as foUows [4] 

r = /ext(r) = Qi(r-m)+Qext^(r) (21) 
where Qi y Qext are arbitrary constants and the function g{r) is defined as 

/ 2m\ 

g{r) = (r - m) log ( 1 —j + 2m (22) 

Looking at the asymptotic condition, we have 

/ 771/ \ ^777^ 

/ext(r ^ oo) ~ Qir 1 - - + Qext Qir (23) 

\ r / r 

and hence the constant Qi must be equal to one (Qi = 1). 

B) The interior problem 

By now introducing expression (1) from the interior metric into the generic equation 
(19), we obtain the foUowing differential equation for the function /(r) 



r^7(r) f"{r) + r 



27(r) - (rp) - 272 (r) 

372 (ro) — 72 (r) 



/'(r)-2/(r) = (24) 



which is at first sight difficult to tackle. Nevertheless, the change of the independent 
variable 

r r e [0 , ro] 

r — >a; = l-72(r) : < (25) 

\xe [0,a;o = 1 -7^(''o)] 

leads to the foUowing equation 
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where primes now denote derivatives with respect to the new variable x and where the 
foUowing definition has been used 



6 = l-37^(ro) 



(27) 



Let US note here what restrictions the upper limit of the parameter /xq imphes: 



This new differential equation (26) has four singular points, aU of which are regulär 



Now, since we are interested in Solutions with good behaviour in the [0, a;o] interval, we 
look for analytical Solutions at the origin as a Probenius series. This procedure leads to a 
Single analytical series whose exponent is Ef = 1/2. However, this must be constructed 
by means of a rather awkward recursive expression. Accordingly, we have chosen to 
reduce equation (26) into the canonical Heun form [5] (see Appendix), which allows us 
to obtain the required Solution by using some properties of that equation. The final 
result for the general and analytical Solution at the origin is the foUowing convergent 
expression in the [0, xq] interval 



Qint being an arbitrary constant and where H{a,q; a, ß,j, ö; z) is the so-called Heun 
series, defined as foUows 




(28) 



{x = b,x = 0,x = 2,x = oo} 



(29) 




(30) 



oo 




(31) 



with the recursive law (n = 0, 1, 2, . . .) 



o(n + 2)(7 + n + l)cn+2 



(n+1) a + ß-S + n + l + a{-f + S + n) + 



n+1 



Cn+l 



+ [77,(0; + ß + n) + aß] Cn = 



= 1 



ci = — 

07 



(32) 
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which is convergent ifF 



(33) 



Note that in our case we have 

2 X 3 

l<|-|<oo , |tI<1 if /^o < 7 (34) 
4 

i.e., the series is convergent everywhere inside the star only if its radius is larger than 
the Schwarzschild radius by a factor 4/3. For example, a typical neutron star of 1.4 solar 
masses and a radius of 2 x 10^ m gives a parameter 

^„ = Hül = = 0.15 < 0.75 (35) 

ro 20 

Thiis, /io < 3/4 is a high upper limit that does not imply restrictions to a realistic star 
model and hence Solution (30) becomes broadly general. 

By writing Solution (30) as a power series in the parameter r/L, we find the foUowing 
expression 



(^ ^ + ^^ 8 + 1% + 456^ 



152 + 4246 + 7296^ + 15756^ 
^ 15120 63 L6 ^ 



(36) 



• Quan's case 



For this case, Solution (30), (36) is strongly divergent and we must return to the original 
differential equation (24), taking the parameter 6 to be zero (37^ (ro) = 1). The resulting 
equation only has three regulär singular points (onc of thcm at infinity) and hence reduces 
to a hypergeometric equation by using the Standard procedure (see [4] for the final result). 

C) Admissible Lichnerowicz coordinates 

We shall prove that it is possible to set the constants Qint and Qext in such a way that 
the components of the metric and their derivatives, written in the new set of coordinates, 
are continuous on the surface of the star (r = ro). In order to do so, we only require the 
function /(r) defining the change of coordinates to be C-^ at that boundary, i.e., 

/ext(r'o) = /int(r-o) 

[/]e = [/']e = ^ <; (38) 

/ext(^o) = //nt(^o) 
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which, in agreement with (21-22), (30), can be translated into the following relations 



with 



ro-m + g{ro) Qext = roh{ro) Qi^t 



/.(r).7/|H,3 + l;4, 1,5,1;^ 



Q 



int 



(39) 



(40) 



The set of equations (39) constitutes an algebraic linear System for the unknown variables 
Qint and QextjWhose Solution is 



Qext — 
Qint = 



(rp - m)[h{ro) + rph'jro)] - roh{ro) 
roh{ro)g'{ro) - [h{ro) + roh'{ro)]g{ro) 

(rp - m)g'{ro) - g{ro) 

roh{ro)g'{ro) - [/i(rp) + rp/i'(rp)]^(ro) 



(41) 



which is well defined, since the denominator of both expressions is trivially different from 
zero. We now givc the firsts terms of the expansion in power series of the //p parameter 
for the values of those constants 



. ^ , 12 4 58 

2mgext = -0 ( 3^ - - ^ 



136 
5005 

1439 



/^o + • ■ ■ 



_ 3 1 2 1 3 

Qint - 1 - ^/^o - Y^/^p - + ;l344qq 



(42) 
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It is Worth noticing that these values have nothing to do with those obtained by Quan for 
the case /Up = 8/9, since our results require an upper limit for the parameter //p < 3/4, 
and hence no comparison can be made. 

Let US now check that this choice of constants does indeed lead to a C-*^ metric, 
written in harmonic polar coordinates, at the boundary r = rp. By making use of the 
change (15-18), we have 

ds'' = gu{r)de+^-^df'' + r\de''+sinHd^^) , [r = r(f)] (43) 

Obviously, all the components of the metric as well as the derivatives of gu, gee and gipip 
are continuous. Regarding the derivative of gff, this turns out to be 



1 



grr 
~l2 



f 



f 



Ii ( [drgrr]^ f'iro) - 2^,,(rp) ) (44) 
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However, taking into account the difFerential equations (20) and (24) for the exterior and 
the interior respectively , we have that 

[n^ = -\%l-\ro)nro) (45) 
and therefore, since the discontinuity of drQrr is given by (13), we finally have that 

[dfQrf]^ = (46) 

which is what we wished to prove. Therefore, the harmonic coordinates obtained with 
the constants (41-42) are admissible coordinates in the sense of Lichnerowicz. 

4. - Interpretation of Quan's exterior constant 
A) Multipole expansion 

As aheady mentioned in the Introduction, the purpose of this Section is to search for 
the physical meaning of the constant Qext associated with the choice of gfobal harmonic 
coordinates given in the previous Section. We first carry out a multipole expansion of 
the exterior Scwarzschild metric in those coordinates, i.e., a power series in the inverse 
of the radial coordinate f. This expansion can be taken, up to Order five, from a previous 
paper [7], where it was introduced for other purposes. Here we write the metric on a 
shghtly different basis that is useful for the foUowing calculations. 

ds^ = T{f) dt^ + A{f) Sij + B{f) mn^ dx'dx^ (47) 

with 

^'=J i^i = Srj^') (48) 

T^Mr) , A^'^ , ß-^-^ [r = r(f)] (49) 

By using these definitions and the results from [7] , we obtain the foUowing expansion 
up to the Order 1/f^ 

„ m Km + 
T = -1 + 2— - 2^ + 2-^ - 2 ZT 

A = l + 2— + — +2— + 2-;^ + --^ + --- (50) 



m 



^ „ —iK + „ —6Km + 



B=-^ + 2 +2 



-9Km'^ + 
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where a new constant K is introduced, which is related with Qext as foUows 



K = -m Qext 



(51) 



B) Linear approximation in harmonic coordinates 

As is well known [11], the Einstein equations can be written in the foUowing way 



(52) 



5"^ = -\/— ^5^*^^ being the contravariant metric density, g = det{gaß) the determinant 
of the metric, and the Landau-Lifshitz energy-momentum pseudotensor. One of 
the procedures used in the past to solve these equations is the so-called perturbative 
postminkowskian algorithm [12]; starting from a System of asymptotically Cartesian 
coordinates, this method consists in looking for a Solution as a formal power series in 
the gravitational constant G: 



^aß = ^«/3 ^ ^ (^n jj 



aß 



n=l 



and taking for the coordinates the harmonicity condition 



(53) 



(54) 



In particular, this procedure has been of great interest for the study of certain aspects 
of gravitational radiation [12]. We are now dealing with a vacuum stationary scenario 
-that is, ^ is independent of time- and hence the equations resulting in the linear 
approximation are as foUows 

A 1)^'^ = 

(55) 

di i)'!" = 

whose general Solution can be written in the foUowing way (see, for instance, [13]) 



^'^ = can + d^'W^ + d^'W^ - V^^dpWP , ( A WA = 0) 



(56) 



where the canonical part 1) is given by 

1 



^ oo ^ +1 



1=0 



m=—l 

+1 



(57) 



1=1 



m=—l 



h^'' = 

V can " 
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and Sp being the Geroch-Hansen static and dynamical multipole moments [14] (up 
to a numerical constant factor), and where the usual notation has been used for spherical 
harmonics. Regarding the "gauge" part of the Solution [terms involving the function 
in (56)], we have the foUowing expressions 

oo ^ +1 
/=0 m=-l 



oo ^ +1 oo +i 



1=1 m=—l 



1=1 



m=—l 



°° 1 , -\ +^ oo ^ +1 



1=0 



m=—l 



1=0 



m=—l 



where C[", HP''^, jp y are new constants (scalar or vectorial) whose meanings have 
not yet been discussed. 

The general Solution of equations (55) with spherical symmetry is obtained from the 
previous expressions (57-58) by restricting ourselves to monopolar terms, i.e, by taking 
into account only Z = in the above series: 



/- 1 



^°L = -4^ (m^M°) 



h^-' =0 

V can " 



^ ^can ~ 



c 



K 



= — ^n^ 



(59) 



In agreement with (56), this expressions lead to the foUowing Solution 



f 1 


= -4^ 


r 


< l^oj 


= d^w"^ 








= _2^ 



(60) 



Finally, the linear order of the metric is obtained by using the known relation 

1 } 



1 i^f^ 



(61) 
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which explicitly turns out to be 





m 




^00 


= 2— 




r 






= djw' 


= 






+ 2 


9jk 


= (^7 



(62) 



K 
,3 J 



K 



where we have taken C = since is a gradient and can therefore be omitted (static 
condition). 

Let US notice that, as it should, metric (62) contains all the linear terms of the 
Schwarzschild metric written in harmonic coordinates (47), (50) shown in the previous 
Section. This detail, which might appear insignificant, turns out to be the relevant point 
for under Standing the meaning of the constant K, as we shall now see. 

C) A spherically Symmetrie "singular" source 

A significant feature of the metric (62) is that it proves to be a Solution of the 
lincarizcd Einstein equations with the foUowing singular energy- momentum tensor (order 
zero) on the right hand side of the equations. 



^T°° = m5(f 



1 



. . 

rplj 



-K 



dkiS{x) 



Thus, the equations can be written as foUows: 



(63) 



A f)^^ = 167rm5(f) 
and we can solve them by using the Poisson Integrals to give 



( 1 



fjOO(f) = -4m 



= +2K 



m 



M=-4 

\x — y\ r 



K 



^dV=-2^ö-+6>W 
\x — y\ 



(64) 



(65) 



which is merely the metric density obtained before (60). 
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Let US now calculate the stress quadrupolar moment of this energy- momentum 
tensor (63) 

0.. -.^ 1 



xPx'^T'U^x= -K 
2 



x^x'^dki5{x) d^x 

(66) 



= K 



We see that the constant K, and therefore Qext, are related to the stress quadrupolar 
moment of the source. In our opinion, this is simply a first argument to understand the 
role played by this kind of constant, which are not present when the Thorne gauge is 
used for the analysis of Standard multipole moments. In the future we hope to develop 
a coherent theory able to completely justify the Interpretation made. 

5. - Global asymptotically Cartesian q-harmonic coor- 
dinates 

In Section 3 we looked for a global System of harmonic coordinates that would be 
asymptotically Cartesian and admissible in the sense on Lichnerowicz. This issue can also 
be solved for a System of the so-called q-harmonic coordinates. These coordinates were 
introduced by Bei [8] to analyze reference frames in General Relativity as congruences 
of time-like curves. Although the foUowing results are (almost all) from other authors, 
we shall show them briefly because they have not been published by those authors. 

Q-harmonic coordinates are only meaningful when associated with the also so-called 
q-harmonic congruences [8]. Since the Schwarzschild metric admits as a q-harmonic 
congruence the Killing time congruence, we restrict ourselves to recalling the definition 
of q-harmonic coordinates for this case. The change of coordinates from the Standard 
ones would be as foUows 

{r, ip} ^{x^}^ {x, y, z} ^ {x^} (67) 
with the q-harmonicity condition 

d,(^y/^g'^djx~^^=0 



(68) 



lim g-ij = 5ij 

r— >oo 



whcrc r = yx^ + + , and g represents the determinant of the spatial metric of the 
quotient space defined by the Killing congruence. 

As in the harmonic case, associated polar coordinates are defined by 

X + iy = fsinÖ e**^ 

(69) 

z = f cos 9 
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f = f(r) , S=e , •fi = v (70) 



r— >oo r 

and this leads to the foUowing differential equation 

d f r'^ df 



dr \ Jgrr dr 



-2^f = (71) 



A) The exterior problem (q— harmonic coordinates) 

In this case the previous differential equation (71) is written as follows 

r (r - 2m)f"{r) + 2(r - ^m)/'(r) - 2/(r) = (72) 
whose general Solution was already obtained by Aguirregabiria [9] 

r ^ fe^t{r) = Ji (r -^m) + J2g{r) (73) 



2 

where Ji and J2 are constants, and the function g{r) is defined by 



9ir) = (^-^"^) yi- V ^^^^ 

with the foUowing asymptotic behavior 

/ 3 m 1 \ 

+ (75) 

Thus, according to (73) we have 

Ji + J2 = 1 (76) 

B) The interior problem (q— harmonic coordinates) 

For this case, the differential equation (71) has the foUowing form 

,2 \ / 3^2 



,2 



1 - ^) /"(O + H 2 - ^ ) /'(r) - 2/(r) = (77) 



whose unique analytical Solution at the origin has already been obtained by Teyssandier 
[10] 

r „ / 1 3 5 



/int(r)-P^F^2'2'2'L2 

(78) 

r f 3 r2 9 5 
L 1ÖL2 "'^ 56L4 ^48L6 ^ 
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P being an arbitrary constant, and where F(a, 6, c, ; x) represents the usual hypergeomet- 
ric function whose first terms of its power expansion are shown for clarity of expression. 

C) Continuity at the boundary S 

As in the harmonic case, the arbitrary constants P and J2 (or Ji) can be fixed by 
imposing that the function /(r) be on the surface of the star, i.e., 



[/]e = [/Ie = : <^ 



^0 - Ml + fi'(^o) J2 = roh{ro) — 



Ji + g'iro) J2 = [/i(ro) + ro/i'(ro)] 



(79) 



with 



1 3 5 
2' 2' 2' L2 



h{r) = F 

These conditions lead to the foUowing expressions for the constants 

(ro - 3m/2) [h{rQ) + rji'iro)] - roh{ro) 



roh{ro)[g'{ro) - 1] - [h{ro) + roh'{ro)\ [g{ro) - (ro - 3m/2)] 



(ro - 3m/2)öf'(ro) - g{ro) 



roh{ro)[g'{ro) - 1] - [h{ro) + roh'{ro)] [g{ro) - (ro - 3r^z/2)] 



(80) 



(81) 



whose expansion in power series of the parameter fxo = 2m /ro = Tq/LP' is given by the 
foUowing expressions 



4m^ J2 = ro 



-8 72 
T+35^° 



29 

1Ö5Ö 



825 



(82) 



P 



1 — Mo H ß(\ H 

80^° 280 



FinaUy, it should be mentioned that Teyssandier [10] has proved, in a similar way to that 
used here with the harmonic coordinates, that these values of the constants provide a 
metric on the surface E, i.e., the q- harmonic coordinates obtained are also admissible 
coordinates in the sense of Lichnerowicz. 



6. - Conclusions 

We have presented a global Solution to Einstein equations for static spherically 
Symmetrie perfect fluid with homogeneous energy density distribution. At the outside 
of matter distribution we have, as it should be, the Schwarzschild space-time, whereas 
inside, our Solution coincides with the well-known Schwarzschild Solution. The novelty 
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of Our approach being that the Solution is written globally in harmonic and q-harmonic 
coordinates. At the boundary surface the Lichnerowicz junction conditions are satisfied, 
which imphes that, both, the metric and its first derivatives are continuous across the 
boundary. This, in turn, imphes that the exterior space-time contains an additional 
paramenter, besides the mass. This parameter being related to the stress quadrupolar 
moment [see (66)]. Obviously this new parameter may be eliminated by means of 
a coordinate transformation. All the above brings out two facts: on one hand, the 
Schwarzschild space-time may be produced by a variety of sources (which is well-known) 
and on the other hand, when the Solution is written globally, using Lichnerowicz junction 
conditions, it contains an additional parameter reflecting properties of the source, other 
than the mass. 

Appendix : 

Heun's equation 

A) The so-called Heun equation is the foUowing second-order linear differential 
equation 

x{x -l){x- a)F"{x) + {{a + ß + l)x'^ 

{AI) 

-[a + ß + l + a{^ + S) - ö]x + a'y}F'{x) + {aßx - q)F{x) = 

where (a, q; a, /3, 7, ö) are numerical constants. This equation has the foUowing singular 
points, all of which are regulär 

{x = 0,x = l,x = a,x = 00} (^2) 

We restrict oursei ves to the point a; = 0, which is the one of interest in the main text, 
and one readily sees that the Frobenius series type Solutions have exponents (0,1 — 7). 

For the zero exponent, the Solution is 

F{x) = H{a,q;a,ß,^,6;x) , (7^0,-1,-2,...) (^3) 

H being the Heun series (31) appearing in the main text, which is convergent if |a| > 1 
and with a radius of convergence < 1. 

For the exponent 1 — 7, the Solution is 

F{x) = x^-^H{a,qi;ai,ßi,^i,ö;x) , (7^1,2,3,...) (^14) 

with 
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51 = 5 + (1 - 7)(a + /? + 1 - 7 - <^ + a<^) 

{Ab) 

ai = a + 1 - j , ßi = ß + 1 - , 7i = 2- 7 

When the previous series are not convergent, i.e., \a\ < 1, we perform the change of 
variable 



ä = - ) {A6) 



a 



and this leads to the following Solution 



Ef = : F{x) = H{a, q; a, ß, 7, ö; x) 

Ef = l-^ : F{x) = x^~^H{ä,qi;ai,ßi,^i,S;x) 



with 



q = aq , S = a + ß+ l — j — S , qi = aqi (^8) 

B) The differential equation (26) appearing in the main text can be reduced to the 
Heun canonical form by means of the following change of function 

/ — >F: f{x) = x^{x-2)\x-h)'^ F{x) {AQ) 

A rather awkward but straightforward calculation shows that there are four sets 
of suitable values for the exponents {k ,1 ,m) and consequently for the corresponding 
Parameters {a , q;a , ß j'y ,S). Since we are looking for an analytical Solution of equation 
(AI) in the neighbourghood of the origin {x = 0), we obtain eight possibilities, provided 
by the two Frobenius exponents (0,1 — 7) of the Heun equation at this point. Taking 
into account that the only correct Frobenius exponent from the starting equation (26) 
is 1/2, we finally have only four possibilities. It can be checked that, naturally, all these 
four possibilities lead to the same analytical Solution at a; = 0, one of them being defined 
by the following set of parameters 

Ä;=- , ^ ~ 2 ' "^ = 
Ef = { {AlO) 
b b 1 , o . 

0=2' ^"^2^2' ' P = l, 7=2' ^=2 

Now, since \a\ < 1, we must resort to formulae (A6), (A8), which provide the following 
values for the parameters 

ä=l , Q = ^+l , « = 4 , ß = l , 1=1 , S=l {AU) 
which construct the Solution shown in the main text (36). 
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